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A METHOD OF SOLVING LINEAR DIFFERENTIAL EQUATIONS. 

SECOND PAPER. 

By P. A. Lambert. 

The object of this paper is to apply to linear partial differential equations 
the method of solution applied to ordinary linear differential equations 
in the paper entitled "A Method of Solving Linear Differential Equations" 
published in the Annals of Mathematics, July, 1910. 

Let the given differential equation be 

The method of solution proposed consists of the following steps: 

(a) Break up the function / into two parts, one of which, /i, equated 
to zero gives a differential equation which may be readily solved, and 
introduce a parameter t as a factor of the second part / 2 , so that the given 
equation, /i + / 2 = 0, is replaced by 

(2) h + tfi = 0. 
(6) Assume that the series 

(3) U = U + Ult + U>£' + u 3 t 3 + • • • , 

where u , u h ik, u 3 , • • • are undetermined functions of x and y, makes 
equation (2) an identity. Substitute the expression (3) in equation (2) 
and determine these functions by solving the differential equations formed 
by equating to zero the coefficients of successive powers of t in the resulting 
identity. 

(c) Substitute these values of u a , u h Ui, u it • • • in (3), and replace t by 
unity. Then see if 

(4) u = u 9 + Ui + Ui + u 3 + • • • 

is convergent and satisfies equation (1). 

This method of solving differential equations will be called the parametric 
method. 

The results obtained by the application of this method are not new. 
The actual solution of an ordinary linear differential equation by the 
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2 P. A. LAMBERT. 

parametric method is simpler in theory and decidedly less laborious than 
by the method which assumes the solution to be 

y = Z) Ajc m+r> 

and requires the determination of the constants A r , m and s. 

The method of solution of linear ordinary differential equations out- 
lined by Schlesinger* is practically identical with a special application of 
the parametric method. In the application of the method outlined by 
Schlesinger, which was established by Caqu6f by studying equations of 
finite differences, the given differential equation, 

/ dy d"y\ 

f { x ' y >dx'---dx") = 
is broken up into 

fi-fc = p(x) or Di(y) - A(y) = p{x), 

where A must contain the derivative of highest order. The solution of 
Di(y) = is called Ua and the solution of the given equation is written 
y = Uo + u, so that 

Di(u + u) = Di(uo + u) + p(x), 

and u must be the principal integral of 

D 1 (u) - A (u) = Fo(x), 

where F<,(x) represents the known function 

p(x) + D 2 (u ). 

If Ui is the principal integral of 

A(w) - F (x) 

and u — Mi + v, v must be the principal integral of 

Di(») = A (f) + ^i(x) 
where 

F,(x) = A(«i). 

Repeated apphcation of this process gives 

y = u + ui + u» + u 3 + ■■-, 



*Handbuch der Theorie der linearen Differentialgleichungen, pp. 370-377. 
tLiouville's Journal, Series II, Vol. 9, p. 185. 
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which is proved to be convergent and a solution of the given differential 
equation. 

In applying the parametric method the given differential equation is 
replaced by 

/i + tf 2 = 0, 

where /i is selected so that (1) the equation /i = can be integrated 
and (2) the resulting series shall be convergent. It is frequently advan- 
tageous to select /i so that it will not contain the highest derivative of the 
given equation. 

Schwarz* discusses the partial differential equation 

dhi dht, 

and Darbouxf discusses the partial differential equation 

d 2 u du , , du , 

■; — r~ = a~ — \- b — + cu, 
dxdy dx dy ' 

where a, b, c, are functions of x and y, by a method of successive approxi- 
mations which is essentially the same as the parametric method. However, 
the use of the parameter, which characterizes my method of solution, 
does not seem to occur in the literature of differential equations. Moreover, 
the parametric method seems better adapted to the actual determination 
of the solution of the partial differential equation than the method of suc- 
cessive approximations as used by Schwarz and Darboux. 

The parametric method will be exemplified by applying it to several 
examples. 

In Example I the method is applied to a first order equation to throw 
the method into prominence. Example II is the second order differential 
equation of fundamental importance in mathematical physics and the 
theory of functions. The solution containing two arbitrary functions is 
found. Examples III, IV and V show how to find particular integrals 
under different conditions 

The differential equations of the first five examples have constant 
coefficients. In Example VI, the differential equation of Euler and Poisson, % 
of importance in mathematical physics and differential geometry, the coeffi- 
cients are functions of x and y. 



*Abhandlungen, Vol. I, pp. 241-265. 

t Theorie generate des surfaces, vol. IV, pp. 353-367. 

% Darboux, Theorie des surfaces, vol. II, p. 54. 
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Example I. 
Replace the equation by 



du , du 

3x + Ty =nu - 



and assume 
There results 



du du, 

nu + -r- t = 

By dx 

U = Uq + Uit + u4 2 + m 3 £ 3 + 



t 2 + 



dUo 


dtk 

dy 




nu 


— nu\ 


— nui 




, du 

+ dx 


, du x 
~ t ~ dx 



= 0. 



Solving the differential equations obtained by equating to zero the 
coefficients of powers of t in this identity, 

Ua = e ny ip{x), 

u 1 =-e""<p'(x)y, 

u. i =e»< P "(x)\ v 



Substituting in the value of u and making t unity, 

u = e"»[<p(x)-v'(x)y + <p"(x)^ J, 

whence by Taylor's series 

u = e ny tp(x — y), 

the integral of the given differential equation containing an arbitrary 

function. 

d 2 u 



Example II. 



dx 



,+« 2 5 = o. 



Replacing this equation by 

d 2 u 



W 



,d 2 u 



and assuming 



TT + Otrr =0, 
dx 2 dy 2 

U = Uo + Uit + Uit 2 + u 3 P + 
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there results 



dx 2 + dx 2 



t + 

+ a 2 



~dx 2 

ey 2 



t 2 + ■ ■ ■ = 0. 



From the series of differential equations obtained by equating to zero 
the coefficients of powers of t in this identity may be found 



uo = (A + Bx)<p(y), 
u^-tf^A^ + B^^'iy), 



Substituting in u and making t unity, 



a'x' 



a 4 x* 



u = A v (y) -A^r f"(y) + A ~4T ^^ + 



.era' 



a*x° 



+ B<p(y)x - B -gy- /(*/) +B-^- v iv (y) . 

Now if <p(y) is made e" in order that <p(y) and all its derivatives shall 
have the same value, and if the arbitrary constant B is replaced by Ca, 
it is readily seen that 

u = Ae" cos (ax)+Ce y sin (ax). 

It is evident that this solution of the differential equation may be 
written in the form 

. „„ ax , „ „„ . ax 

u = Ae v cos h Ce " sin — . 

n n 

If it is desired that u shall contain all the successive derivatives of <p(y) 
the values of u 0) u\, u%, Uz, • •• may be written 

«o = <p(y) + Bx<p'(y) 



2-2 



a l x 



2„S 



erx 



mi = -<P H (y)-2\-B v '"(y)-^ 
U2 = <f> lv (y) -jf + B<p v (y) -jjy 



a'x' 
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Whence 

« = <p(y) + <p'(y)Bx - < P "(y)' z ^ 

n 2 r 3 ,,4-4 

-/"(2/)5^+^(2/)- 4 y+---- 

The terms of n which have been written are the first five terms of the 
expansion of a function by Taylor 's series provided B = ai or B - — ai, 
where i — V — f. 

The terms of u which have been written now suggest that 

u = <pi(y + tax) + <p 2 (y - iax), 

which on trial is found to be correct. 

Of this solution containing two arbitrary functions the former solution 
containing two arbitrary constants is a special case. 



Example 


III. 


dx 2 dy ~ 


2 d ' y + 


a 3 w 
dy* ~ 


1 

x 2 


Replacin 


Lg this 


equation by 












d s u 1 


.(*-** 


d 3 U 


\t 



dx*dy x 2 \ dxdy 2 dy 3 J 

and applying the method 

Mo = <pi(x) + (A + Bx)<f> 2 (y) + (Ax + B!x)^(y) - y log x. 

By actual trial it is found that the following four values of u 0) parts of 
this general expression for u , 

u a = <pi(x), 

«o = <P2(y) + Bix<p'. 2 {y), 

u a = Bx<p z (y), 

w = - y log x, 

lead to solutions of the differential equation. 

The solution of the equation, the sum of the solutions corresponding to 
these four values of u , is 

u = ^x(x) + i&ix + y) + x<p 3 {x + y) - y log x. 
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Example IV. - 2 - -, - 3 - + 3 j- y - ^ + a*. 

The complementary integral is the solution of the equation 

dx 2 ~ dy 2 ~ 6 dx + 6 dy~ °' 
Replacing this equation by 



d 

dx 



u du (d^u „du\ 

: 2 6 dx \dy 2 ~ 6 dy) t ~ [J ' 



we find u = e Zx <p(y), which suggests that e 3x is a factor of the comple- 
mentary integral. 

Transforming the equation by the relation u = e 3x -v, a partial differ- 
ential equation is formed of which 

v = <p(x - y) 

is a solution. It follows that 

u = e 3x <p(x — y) 

is a complementary integral. 

If the differential equation which determines the complementary 
integral is replaced by 



(eju_ du\ (dhi_ au\_ 



it is found in like manner that 

u — e 3 "<p(x — y) 

is a complementary integral. 

These two complementary integrals are not independent. Replacing 
the differential equation by 

( d 2 U _ d 2 u\ _ (du du\ 

\dx 2 dy 2 ) 6 \dx ay)~ ' 
a third and independent complementary integral is found 

u = <p(x + y). 
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The particular integral corresponding to the term xy is found by applying 
the method to 

b 2 u 
dx 2 



du (b 2 u du\. 



To find the particular integral corresponding to the term e i+2y transform 
the equation 



d 2 u d 2 u „du . „ du 
by the relation 



dx 2 by 2 dx ^ ° by ~ 



u = ve' +2y . 
The equation in v is 



d 2 v _ d 2 v _ dv _ dv 
dx 2 ~ dy 2 dx dy ~ ' 



A special solution of this equation is v = — y, and the corresponding 
particular integral of the given equation is 

u = — y-e x+2y . 

The solution of the given differential equation is the sum of the two 
independent complementary integrals and the two parts of the particular 
integral. 

Example V. To find a particular integral of 

d-u , d 2 u „d 2 u j . , . , 

ax- 2 + 8x^- Q W t==XSm(x + y) - 

Replace the right-hand member by x 2 e" +iy and transform the resulting 
equation by the relation 

u = e iz+ty -v. 

The particular integral of the equation in v can be found by the method 
of example IV. 

The particular integral of the given equation is the coefficient of i in 
the corresponding value of u. 

So far the parametric method has been applied to linear partial differ- 
ential equations with constant coefficients. The following example shows 
that the method may be applied with advantage to equations with variable 
coefficients. 

„„ d 2 u a du , b du n 

Example VI. 7— H r- = 0. 

*«""v * dzdy x—ydx x — y dy 
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du , dU . . d"U , 

a~ o-r- = (x — y) r— — t 

dx dy v *" ax ay 



Replacing this equation by 

and assuming that 

U = U + Uit + Uit 2 + u 3 t 3 + • • • , 

it follows at once that 

Mo = <p(ay + bx). 

Solutions of the given equations which are homogeneous polynomials 
are found by giving u the successive values 

u = ay + bx, 

Mo = (ay + bx)-, 

Mo = (ay + bx) 3 , 



The given equation may be written 
/ a 2 M 



x- 



. ,du\ ( d 2 u du\ 

+ b Ty)-{ y ^xTy + a Tx)=°- 



dxdy 

A solution of the first part of this equation placed equal to zero is 

u = x - Vi(w) ; 

a solution of the second part placed equal to zero is 

m = y~ a <p2(x). 

Hence a solution of the given equation is 

m = x- b y~ a , 

from which we may get the more general solution 

u = (x + m)-\y + m)- a , 

where m is an arbitrary constant. 

The given differential equation may be replaced by either of the equations 

d 2 M a du b 3m . _ n 



dxdy x — y dx x — ydy 
d 2 u . b du a du _ 



dxdy x — ydy x — y dx 



10 P. A. LAMBERT. 

From the first equation a special value of m is 



tin = 



I - a 



The use of this value of u suggests that (x — y) 1-a is a factor of the solution 
of the given equation. In like manner from the second equation 

-Or- y y- b 
uo = — jrrr - , 

which suggests that (x — y) l ~ h is a factor of the solution of the given 
equation. 

From these suggestions it is inferred that (x — y) 1 ~ a ~ b is a factor of 
the solution of the given equation. This inference is correct, for the relation 

u = (x — y) l - a ~ h -v 

transforms the given equation into 

d-v 1 — b dv l_^_a dv _ 



dxdy x — y dx x — y dy 

By the preceding paragraph a solution of this equation is 

v = (x + my-^y + n) h -\ 

The corresponding solution of the given equation is 

u = (x + m)°~\y + mf-^x - j/) 1 ""- 6 . 
Lehigh University, February, 1911. 



